Efficiency of membrane fusion inhibitors on different hemagglutinin subtypes: insight from a molecular dynamics simulation perspective.

The challenge in vaccine development, along with drug resistance issues, has encouraged the search for new anti-influenza drugs targeting different viral proteins. Hemagglutinin (HA) glycoprotein, crucial in the viral replication cycle, has emerged as a promising therapeutic target. CBS1117 and JNJ4796 were reported to exhibit similar potencies against infectious group 1 influenza, which included H1 and H5 HAs; however, their potencies were significantly reduced against group 2 HA. This study aims to explore the molecular binding mechanisms and group specificity of these fusion inhibitors against both group 1 (H5) and group 2 (H3) HA influenza viruses using molecular dynamics simulations. CBS1117 and JNJ4796 exhibit stronger interactions with key residues within the H5 HA binding pocket compared to H3-ligand complexes. Hydrogen bonding and hydrophobic interactions involving residues, such as H381, Q401, T3251 (H5-CBS1117), T3181 (H5-JNJ4796), W212, I452, V482, and V522 predominantly contribute to stabilizing H5-ligand systems. In contrast, these interactions are notably weakened in H3-inhibitor complexes. Predicted protein-ligand binding free energies align with experimental data, indicating CBS1117 and JNJ4796's preference for heterosubtypic group 1 HA binding. Understanding the detailed atomistic mechanisms behind the varying potencies of these inhibitors against the two HA groups can significantly contribute to the development and optimization of effective HA fusion inhibitors. To accomplish this, the knowledge of the transition of HA from its pre- to post-fusion states, the molecular size of ligands, and their potential binding regions, could be carefully considered.Communicated by Ramaswamy H. Sarma.

Insights into binding molecular mechanism of hemagglutinin H3N2 of influenza virus complexed with arbidol and its derivative: A molecular dynamics simulation perspective.

Recently, the H3N2 influenza outbreak has caused serious global public health concern for future control of the next influenza pandemic. Since using current anti-influenza drugs targeting neuraminidase (oseltamivir and zanamivir) and the proton M2 channel (amantadine and rimantadine) leads to drug resistance, it is essential to seek new anti-viral agents that act on additional viral targets. Hemagglutinin (HA), a glycoprotein embedded in the viral surface and playing a critical role in influenza the viral replication cycle has become an attractive target. This work investigates the molecular binding mechanism of HA H3N2 of influenza virus complexed with the fusion inhibitor, arbidol and its derivative (der-arbidol), by means of molecular dynamics simulation. The result showed that the arbidol derivative could form many and strong hydrogen bonds with the HA surrounding amino acids comprising GLU1032(1), LYS3071(1) and LYS3102(1) while arbidol makes this type of interaction with only LYS582(1). The introduction of hydroxyl group at the meta-position of the thiophenol ring was detected to replace the nearby water molecule, thus allowing the direct hydrogen bond formation between der-arbidol and GLU1032(1) of HA residue. Furthermore, the salt bridge networks established among residues GLU572(1)···ARG542(1)···GLU972(2) were considerably more stable in HA-Der-arbidol than that found in HA-Arbidol. The predicted protein-ligand binding free energies were in agreement with experimental data indicating that der-arbidol exhibits higher inhibitory potency against HA H3N2 of influenza virus. Detailed information could be useful for further designing and optimizing HA fusion inhibitors with improved efficiency.

Exploring of paritaprevir and glecaprevir resistance due to A156T mutation of HCV NS3/4A protease: molecular dynamics simulation study.

Hepatitis C virus (HCV) NS3/4A serine protease is a promising drug target for the discovery of anti-HCV drugs. However, its amino acid mutations, particularly A156T, commonly lead to rapid emergence of drug resistance. Paritaprevir and glecaprevir, the newly FDA-approved HCV drugs, exhibit distinct resistance profiles against the A156T mutation of HCV NS3/4A serine protease. To illustrate their different molecular resistance mechanisms, molecular dynamics simulations and binding free energy calculations were carried out on the two compounds complexed with both wild-type (WT) and A156T variants of HCV NS3/4A protease. QM/MM-GBSA-based binding free energy calculations revealed that the binding affinities of paritaprevir and glecaprevir towards A156T NS3/4A were significantly reduced by ∼4 kcal/mol with respect to their WT complexes, which were in line with the experimental resistance folds. Moreover, the relatively weak intermolecular interactions with amino acids such as H57, R155, and T156 of NS3 protein, the steric effect and the destabilized protein binding surface, which is caused by the loss of salt bridge between R123 and D168, are the main contributions for the higher fold-loss in potency of glecaprevir due to A156T mutation. An insight into the difference of molecular mechanism of drug resistance against the A156T substitution among the two classes of serine protease inhibitors could be useful for further optimization of new generation HCV NS3/4A inhibitors with enhanced inhibitory potency.Communicated by Ramaswamy H. Sarma.

Equivalence between inverted regions of the energy gap law and inverted regions of donor-acceptor distances in photoinduced electron transfer processes in flavoproteins.

In the present work, we discuss about the relationship between the energy gap law and extended Dutton law in flavoproteins. The extend Dutton law is defined herein as the dependence of logarithmic rates (ln Rate) of photoinduced electron transfer (ET) from aromatic amino acids to excited isoalloxazine (Iso*) on donor-acceptor distances (Rcs). Both functions of ln Rate vs. negative values of the standard free energy gap and ln Rate vs. Rc display a parabolic behavior, when the ET rates are ultrafast. The negative values of the standard free energy gap at peaks of ln Rate [X m(ES)] were obtained for FMN-binding protein, wild-type pyranose 2-oxidase, T169S (Thr169 is replaced by Ser) pyranose 2-oxidase, and medium-chain acyl-CoA dehydrogenase. The values of Rc at peaks of ln Rate [X m(Rc)] were also obtained for these flavoproteins. The negative values of the standard free energy gap decreased with approximate linear functions of Rc. The negative values of standard free energy gap [X m(ESRc)] at Rc = X m(Rc) were evaluated using the linear functions of the negative standard free energy gap with Rc. The values of X m(ESRc) were mostly in very good agreement with the values of X m(ES). This implies that the energy gap law and the extend Dutton law are equivalent. X m(ES) values in ET donors displaying the linear extend Dutton law with Rc were obtained by energy gap law, and then X m(Rc) values were evaluated with the negative standard free energy gap. Thus, the obtained X m(Rc) values were much smaller than the Rc range obtained by the method of molecular dynamics simulation. This suggests that ET processes with linear profiles of the extend Dutton law could be parabolic when Rc becomes much shorter than the Rc range obtained by the method of molecular dynamics simulation.

In Silico ADME, Metabolism Prediction and Hydrolysis Study of Melatonin Derivatives.

Melatonin (MLT) is a well-known pineal hormone possessed with remarkable biological activities. However, its low oral bioavailability and high first-pass metabolism rate are important pharmacokinetics problems. Therefore, 5 MLT derivatives (1-5) were designed and synthesised in our group to solve these problems. In this work, in silico analysis of all synthetic derivatives for pharmacokinetic and drug-likeness parameters were predicted by SwissADME software. The results revealed that all derivatives (1-5) met the requirements for ideal oral bioavailability and CNS drugs. The molecular docking showed that the acetyl-MLT derivative (1) and the un-substitution at N1-position derivative 5 would be substrates of CYP1A2, while the lipophilic substituted N1-position derivatives 2-4 could not be metabolised by CYP1A2. Moreover, all N-amide derivatives (1-4) were hydrolysed and released less than 2.33% MLT after 4-hour incubation in 80% human plasma. It seemed that these derivatives preferred to behave like drugs rather than prodrugs of MLT. These findings confirmed that the addition of bulky groups at the N1-position of the MLT core could prolong the half-life, increase drug absorption and penetrate the blood brain barrier into the CNS.

Understanding of the drug resistance mechanism of hepatitis C virus NS3/4A to paritaprevir due to D168N/Y mutations: A molecular dynamics simulation perspective.

Hepatitis C virus (HCV) NS3/4A protease is an attractive target for the development of antiviral therapy. However, the evolution of antiviral drug resistance is a major problem for treatment of HCV infected patients. Understanding of drug-resistance mechanisms at molecular level is therefore very important for the guidance of further design of antiviral drugs with high efficiency and specificity. Paritaprevir is a potent inhibitor against HCV NS3/4A protease genotype 1a. Unfortunately, this compound is highly susceptible to the substitution at D168 in the protease. In this work, molecular dynamics simulations of paritaprevir complexed with wild-type (WT) and two mutated strains (D168 N and D168Y) were carried out. Due to such mutations, the inhibitor-protein hydrogen bonding between them was weakened and the salt-bridge network among residues R123, R155 and D168 responsible for inhibitor binding was disrupted. Moreover, the per-residue free energy decomposition suggested that the main contributions from key residues such as Q80, V132, K136, G137 and R155 were lost in the D168 N/Y mutations. These lead to a lower binding affinity of paritaprevir for D168 N/Y variants of the HCV NS3/4A protease, consistent with the experimental data. This detailed information could be useful for further design of high potency anti-HCV NS3/4A inhibitors.

Physical quantity of residue electrostatic energy in flavin mononucleotide binding protein dimer.

The electrostatic (ES) energy of each residue was for the first time quantitatively evaluated in a flavin mononucleotide binding protein (FBP). A residue electrostatic energy (RES) was obtained as the sum of the ES energies between atoms in each residue and all other atoms in the FBP dimer using atomic coordinates obtained by a molecular dynamics (MD) simulation. ES is one of the most important energies among the interaction energies in a protein. It is determined from the RES, the residues which mainly contribute to stabilize the structure of each subunit, and the binding energy between two subunits can be estimated. The RES of all residues in subunit A (Sub A) and subunit B (Sub B) were attractive forces, even though the residues contain net negative or positive charges. This reveals that the ES energies of any of the residues can contribute to stabilize the protein structure. The total binding ES energy over all residues among the subunits was distributed between -0.2 to -1.2 eV (mean = -0.67 eV) from the MD simulation time.

Dynamics of the protein structure of T169S pyranose 2-oxidase in solution: Molecular dynamics simulation.

Pyranose 2-oxidase (P2O) from Trametes multicolor contains FAD as cofactor, and forms a tetramer. The protein structure of a mutated P2O, T169S (Thr169 is replaced by Ser), in solution was studied by means of molecular dynamics simulation and analyses of photoinduced electron transfer (ET) from Trp168 to excited isoalloxazine (Iso*), and was compared with wild type (WT) P2O. Hydrogen bonding between Iso and nearby amino acids was very similar as between T169S and WT protein. Distances between Iso and Tyr456 were extremely heterogeneous among the subunits, 1.7 (1.5 in WT) in subunit A (Sub A), 0.97 (2.2 in WT) in Sub B, 1.3 (2.1 in WT) in Sub C, 1.3 nm (2.0 in WT) in Sub D. Mean values of root of mean square fluctuation over all residues were greater by four times than those in WT. This suggests that the protein structure of T169S is much more flexible than that of WT. Electrostatic (ES) energies between Iso anion in one subunit and ionic groups in the entire protein were evaluated. It was found that more than 50% of the total ES energy in each subunit is contributed from other subunits. Reported fluorescence decays were analyzed by a method as WT, previously reported. Electron affinities of Iso* in T169S were appreciably higher than those in WT. Static dielectric constants near Iso and Trp168 were also quite higher in T169S than those in WT.

Conformational difference between two subunits in flavin mononucleotide binding protein dimers from Desulfovibrio vulgaris (MF): molecular dynamics simulation.

The structural and dynamical properties of five FMN binding protein (FBP) dimers, WT (wild type), E13K (Glu13 replaced by Lys), E13R (Glu13 replaced by Arg), E13T (Glu13 replaced by Thr) and E13Q (Glu13 replaced by Gln), were investigated using a method of molecular dynamics simulation (MDS). In crystal structures, subunit A (Sub A) and subunit B (Sub B) were almost completely equivalent in all of the five FBP dimers. However, the predicted MDS structures of the two subunits were not equivalent in solution, revealed by the distances and inter-planar angles between isoalloxazine (Iso) and aromatic amino acids (Trp32, Tyr35 and Trp106) as well as the hydrogen bonding pairs between Iso and nearby amino acids. Residue root of mean square fluctuations (RMSF) also displayed considerable differences between Sub A and Sub B and in the five FBP dimers. The dynamics of the whole protein structures were examined with the distance (RNN) between the peptide N atom of the N terminal (Met1) and the peptide N atom of the C terminal (Leu122). Water molecules were rarely accessible to Iso in all FBP dimers which are in contrast with other flavoenzymes.

Classification of the mechanisms of photoinduced electron transfer from aromatic amino acids to the excited flavins in flavoproteins.

In many flavoproteins photoinduced electron transfer (ET) efficiently takes place from aromatic amino acids such as tryptophan or tyrosine to the excited isoalloxazine, so that the fluorescence lifetimes of isoalloxazine in some flavoproteins become ultrashort. The mechanism of ET in the flavoproteins was classified into four classes from the relationship between logarithmic ET rates (ln Rate) and the donor-acceptor distances (Rc), using reported data. The physical quantity, GT, is defined as the sum of solvent reorganization energy, electrostatic energy between a donor cation and an Iso anion, the standard free energy gap between the photoproducts and reactants, and net electrostatic energy between the photoproducts and other ionic groups in the flavoproteins (NetES). When GT fluctuates around zero with Rc, the ET rate becomes fastest (faster than 1 ps(-1)) in Kakitani and Mataga rates. In the ultrafast ET processes, the ln Rate becomes a parabolic function (category 1) of Rc as in FMN binding proteins and pyranose 2-oxidase at the shorter emission wavelengths, when NetES is negligible compared to the other quantities in the GT function. In the ultrafast ET processes, the ln Rate does not display any clear function of Rc (category 2) when NetES is dominant in the GT function, because of no direct relation between NetES and Rc. ET in flavodoxin from Helicobacter pylori may be classified into category 2. When GT linearly varies with Rc around a certain positive value, the ET rates become much slower (<1 ps(-1)). In this case the ln Rate linearly decreases with Rc (category 3), as Tyr224 in d-amino acid oxidase dimers. It is also conceivable that the ln Rate decreases with much scattered function of Rc (category 4), when NetES is dominant in the GT function, as Tyr314 in d-amino acid oxidase dimers. In ET processes of category 1, ET rates decrease as Rc becomes shorter than the distance at the maximum values of ln Rates, where GT is negative. Conditions and physical meanings were discussed for the GT-negative region.

Non-equivalent conformations of D-amino acid oxidase dimer from porcine kidney between the two subunits. Molecular dynamics simulation and photoinduced electron transfer.

The structural difference between two subunits of D-amino acid oxidase dimer from porcine kidney was studied by molecular dynamics simulation (MDS) and rate of photoinduced electron transfer (ET) from aromatic amino acids as tyrosines (Tyr) and tryptophanes (Trp) to the excited isoalloxazine (Iso*). The donor-acceptor distances (Rc) between isoalloxazine (Iso) and the donors were shortest in Tyr224 (0.74 nm) in Sub A at 10 °C (Sub A10), in Tyr224 (0.79 nm) in Sub B at 10 °C (Sub B10), in Tyr228 (0.85 nm) in Sub A at 30 °C (Sub A30), and in Tyr224 (0.72 nm) in Sub B at 30 °C (Sub B30). The Rcs were mostly shorter in the dimer than those in the monomer. Hydrogen bonding (H-bond) pairs between Iso and surrounding amino acids varied with the subunit and temperature. O2 of the Iso ring formed an H-bond exclusively with Thr317OG1 (side chain) in both Sub A10 and Sub A30, while it formed with Gly315N (peptide), Leu316N and Thr317N in Sub B10 and Sub B30. N3H of Iso formed an H-bond with Leu51O (peptide) in Sub A10 and Sub A30, but not in Sub B10 and Sub B30. Electron affinity of Iso* was appreciably lower in Sub A10 compared to Sub B10, while it was opposite at 30 °C. ET rate to Iso* was fastest from Tyr224 in Sub A10, while it was fastest from Tyr314 in Sub B10. The ET rate was fastest from Tyr314 in Sub A30, while it was fastest from Tyr224 in Sub B30. The greater ET rates in the dimer as compared to those in the monomer were elucidated with shorter Rc in the dimer as compared to the monomer. The static dielectric constants inside the subunits and the static dielectric constant between Iso and Tyr224 or Tyr228 were not different appreciably. A few water molecules and sometimes an amino acid were located between Iso and Tyr224, which may be the reason why the dielectric constant of the entire subunits did not differ from that between Iso and Tyr224.

Photoinduced electron transfer modeling to simulate flavoprotein fluorescence decay.

A method of analysis is described on the photoinduced electron transfer (PET) from aromatic amino acids as tryptophans (Trp) and tyrosines (Tyr) to the excited isoalloxazine (Iso*) in FMN-binding proteins (FBP) from Desulfovibrio vulgaris (strain, Miyazaki F). Time-dependent geometrical factors as the donor-acceptor distances are determined by means of a molecular dynamics simulation (MDS) of the proteins. Fluorescence decays of the single mutated isoforms of FBP are used as experimental data. The electrostatic (ES) energy between the photoproducts and ionic groups in the proteins is introduced into the Kakitani and Mataga (KM) model, which is modeled for an electron transfer process in solution. The PET parameters contained in the KM rate are determined by means of a nonlinear least square method, according to the Marquardt algorithm. The agreement between the observed and calculated decays is quite good, but not optimal. Characteristics on PET in flavoproteins, obtained by the present method, are described. Possible improvements of the method are discussed.

Binding mode prediction of biologically active compounds from plant Salvia Miltiorrhiza as integrase inhibitor.

Integrase (IN), an essential enzyme for HIV-1 replication, has been targeted in antiretroviral drug therapy. The emergence of HIV-1 variants clinically resistant to antiretroviral agents has lead to the development of alternative IN inhibitors. In the present work, binding modes of a high potent IN inhibitor, M522 and M532, within the catalytic binding site of wild type (WT) IN were determined using molecular docking calculation. Both M522 and M532 displayed similar modes of binding within the IN putative binding pocket and exhibited favorable interactions with the catalytic Mg(2+) ions, the nearby amino acids and viral DNA through metal-ligand chelation, hydrogen bonding and π-π stacking interactions. Furthermore, the modes of action of these two compounds against the mutated Y212R, N224H and S217H PFV IN were also predicted. Although the replacement of amino acid could somehow disturb inhibitor binding mode, almost key interactions which detected in the WT complexes were fairly conserved. Detailed information could highlight the application of M522 and M532 as candidate IN inhibitors for drug development against drug resistant strains.

Effects of the protonation state of the catalytic residues and ligands upon binding and recognition in targeted proteins of HIV-1 and influenza viruses.

The determination of the protonation state of the functional groups of ligands, and the amino acid residues with electrically charged side chains (His, Lys, Arg, Asp and Glu) or the nucleotide bases of the nucleic acids that they interact with, is important for ligand binding and recognition, the enzyme activity and reaction mechanism, and protein folding/unfolding and stability. Herein, the effects of different protonation state assignments of the small substrate and inhibitors and the critical residues on the reverse transcriptase and protease of human immunodeficiency virus type 1 (HIV-1) and the M2 proton channel of influenza A virus are reviewed. Theoretical studies on these topics are summarized and compared with the experimental data.

Binding pattern of the long acting neuraminidase inhibitor laninamivir towards influenza A subtypes H5N1 and pandemic H1N1.

Influenza A H5N1 and pH1N1 viruses have broadly emerged and become widespread in various countries around the world. Oseltamivir, the most commonly used antiviral drug against the seasonal and pandemic influenza viruses, is targeted at the viral neuraminidase (NA), but some isolates of this virus have become highly resistant to this drug. The novel long-acting drug, laninamivir, was recently developed to inhibit influenza A and B viruses of either the wild-type (WT) or the oseltamivir resistant mutant of NA. To understand the high efficiency of laninamivir, all-atom molecular dynamics simulations were performed on the WT and H274Y mutant of H5N1 and pH1N1 NAs with laninamivir bound. As a result, the novel drug was found to directly interact with 11 binding residues mainly through salt bridge and hydrogen bond formation (as also seen by electrostatic contribution). These are comprised of 7 of the catalytic residues (R118, D151, R152, R224, E276, R292 and R371), and 4 of the framework residues (E119, W178, E227 and E277). Laninamivir showed a similar binding pattern to all four NAs, but strong hydrogen bonding interactions were only found in the WT strain, with a slightly lowered contribution at some drug contact residues being observed in the H274Y mutation. This is in good agreement with the experimental data that the H274Y mutant has a small increase (1.3-7.5-fold, which was not statistically significant) in the IC50 value of laninamivir.

Density functional investigation of hydrogen gas adsorption on Fe-doped pristine and Stone-Wales defected single-walled carbon nanotubes.

The adsorptions of hydrogen molecule of the Fe - doped pristine and Stone - Wales defected armchair (5,5) single - walled carbon nanotubes (SWCNTs) compared with the pristine SWCNT were investigated by using the density functional theory at the B3LYP/LanL2DZ level. The doping of Fe atom into SWCNTs occurring via an exothermic process was found. The adsorptions of hydrogen molecule on the Fe - doped structures of either perfect or SW defected SWCNTs are stronger than on their corresponding undoped structures. The structural and electronic properties of the pristine and SW defected SWCNTs, their Fe - doped structures and their hydrogen molecule adsorptions are reported.

Structural basis for the temperature-induced transition of D-amino acid oxidase from pig kidney revealed by molecular dynamic simulation and photo-induced electron transfer.

The structural basis for the temperature-induced transition in the D-amino acid oxidase (DAAO) monomer from pig kidney was studied by means of molecular dynamic simulations (MDS). The center to center (Rc) distances between the isoalloxazine ring (Iso) and all aromatic amino acids (Trp and Tyr) were calculated at 10 °C and 30 °C. Rc was shortest in Tyr224 (0.82 and 0.88 nm at 10 and 30 °C, respectively), and then in Tyr228. Hydrogen bonding (H-bond) formed between the Iso N1 and Gly315 N (peptide), between the Iso N3H and Leu51 O (peptide) and between the Iso N5 and Ala49 N (peptide) at 10 °C, whilst no H-bond was formed at the Iso N1 and Iso N3H at 30 °C. The H-bond of Iso O4 with Leu51 N (peptide) at 10 °C switched to that with Ala49 N (peptide) at 30 °C. The reported fluorescence lifetimes (228 and 182 ps at 10 and 30 °C, respectively) of DAAO were analyzed with Kakitani and Mataga (KM) ET theory. The calculated fluorescence lifetimes displayed an excellent agreement with the observed lifetimes. The ET rate was fastest from Tyr224 to the excited Iso (Iso*) at 10 °C and from Tyr314 at 30 °C, despite the fact that the Rc was shortest between Iso and Tyr224 at both temperatures. This was explained by the electrostatic energy in the protein. The differences in the observed fluorescence lifetimes at 10 and 30 °C were ascribed to the differences in electron affinity of the Iso* at both temperatures, in which the free energies of the electron affinity of Iso* at 10 and 30 °C were -8.69 eV and -8.51 eV respectively. The other physical quantities related to ET did not differ appreciably at both temperatures. The electron affinities at both temperatures were calculated with a semi-empirical molecular orbital method (MO) of PM6. Mean calculated electron affinities over 100 snapshots with 0.1 ps intervals were -7.69 eV at 10 °C and -7.59 eV at 30 °C. The difference in the calculated electron affinities, -0.11 eV, was close to the observed difference in the free energies, -0.18 eV. The present quantitative analysis predicts that the highest ET rate can occur from a donor with longer donor-acceptor distance, which was explained by differences in electrostatic energy.

Computational studies of influenza A virus at three important targets: hemagglutinin, neuraminidase and M2 protein.

While the seasonal influenza viruses spreading around the world cause the annual epidemics, the recent outbreaks of influenza A virus subtype H5N1 and pandemic H1N1 have raised a global human health concerns. In this review, the applicability of computational techniques focused on three important targets in the viral life cycle: hemagglutinin, neuraminidase and M2 proton channel are summarized. Protein mechanism of action, substrate binding specificity and drug resistance, ligand-target interactions of substrate/inhibitor binding to these three proteins either wild-type or mutant strains are discussed and compared. Advances on the novel anti-influenza agents designed specifically to combat the avian H5N1 and pandemic H1N1 viruses are introduced. A better understanding of molecular inhibition and source of drug resistance as well as a set of newly designed compounds is greatly useful as a rotational guide for synthetic and medicinal chemists to develop a new generation of anti-influenza drugs.

Simultaneous analyses of photoinduced electron transfer in the wild type and four single substitution isomers of the FMN binding protein from Desulfovibrio vulgaris, Miyazaki F.

The mechanism of photoinduced electron transfer (PET) from the aromatic amino acids (Trp32, Tyr35 and Trp106) to the excited flavin mononucleotide (FMN) in the wild type (WT) and four single amino acid substitution isomers (E13T, E13Q, W32A and W32Y) of FMN binding protein (FBP) from the Desulfovibrio vulgaris (Miyazaki F) were simultaneously analyzed (Method A) with the Marcus-Hush (MH) theory and Kakitani-Mataga (KM) theory using ultrafast fluorescence dynamics of these proteins. In addition, the PET mechanism of the WT, E13T and E13Q FBP systems (Method B) were also analyzed with both MH and KM theories. The KM theory could describe all of the experimental fluorescence decays better than the MH theory by both Methods A and B. The PET rates were found to largely depend on the electrostatic energies between photo-products, isoalloxazine (Iso) anion and the PET donor cations, and the other ionic groups, and hence on static dielectric constants. The dielectric constant (ε(0)(DA)) around the PET donors and acceptor was separately determined from those (ε(0)(j), j = WT, E13T, E13Q, W32Y and W32A) in the domain between the Iso anion or the donor cations and the other ionic groups in the proteins. The values of ε(0)(DA) were always lower than those of ε(0)(j), which is reasonable because no amino acid exists between the PET donors and acceptor in all systems. The values of the dielectric constants ε(0)(j) (j = WT, E13T and E13Q) were similar to those obtained previously from the analysis of the crystal structures and the average lifetimes of these FBP proteins. Energy gap law in the FBP systems was examined. An excellent parabolic function of the logarithms of the PET rates was obtained against the total free energy gap. The PET in these FBP isomers mostly took place in the so-called normal region, and partly in the inverted region.

Evaluating how rimantadines control the proton gating of the influenza A M2-proton port via allosteric binding outside of the M2-channel: MD simulations.

In order to understand how rimantadine (RMT) inhibits the proton conductance in the influenza A M2 channel via the recently proposed "allosteric mechanism", molecular dynamics simulations were applied to the M2-tetrameric protein with four RMTs bound outside the channel at the three protonation states: the 0H-closed, 1H-intermediate and 3H-open situations. In the 0H-closed state, a narrow channel with the RMT-Asp44-Trp41 H-bond network was formed, therefore the water penetration through the channel was completely blocked. The Trp41-Asp44 interaction was absent in the 1H-intermediate state, whilst the binding of RMT to Asp44 remained, which resulted in a weakened helix-helix packing, therefore the channel was partially prevented. In the 3H-open state it was found that the electrostatic repulsion from the three charged His37 residues allowed the Trp41 gate to open, permitting water to penetrate through the channel. This agreed well with the potential of the means force which is in the following order: 0H > 1H > 3H.